ON THE ABSOLUTE CONTINUITY OF MULTIDIMENSIONAL 
ORNSTEIN-UHLENBECK PROCESSES 



THOMAS SIMON 

Abstract. Let AT be a n-dimensional Ornstcin-Uhlenbeck process, solution of the S.D.E. 

dXt = AXtdt + dBt 

where A is a real n x n matrix and B a Levy process without Gaussian part. We show 
that when A is non-singular, the law of Xi is absolutely continuous in R" if and only if the 
jumping measure of B fulfils a certain geometric condition with respect to A, which we call 
the exhaustion property. This optimal criterion is much weaker than for the background 
driving Levy process B, which might be very singular and sometimes even have a one- 
dimensional discrete jumping measure. It also solves a difficult problem for a certain class 
of multivariate Non-Gaussian infinitely divisible distributions. 



1. Introduction and statement of the result 

If Z is a real Levy process without Gaussian part, finding a necessary and sufficient 
condition on its jumping measure u for the absolute continuity of Zt at some given t > 0, is a 
hard problem for which no sensible conjecture has been formulated as yet. One of the main 
difficulties for this formulation stems from the time-dependency of the absolute continuity 
property: if u is infinite and has discrete support, then there are some situations where e.g. 
Zi is singular and Z2 is absolutely continous. We refer to [12] and Chapter 27 in [TT] for 
more on this topic as well as further references. When Z is multidimensional, the problem 
becomes increasingly complicated and some partial results had been given in pS] , involving 
conditions of geometrical nature on u. 

On the other hand, the problem of absolute continuity may well become simpler, and yield 
weaker conditions, when considering certain functionals of Z. In the real case for example, 
one can show that 

(1.1) / Ztdtisa.c. <^==^ z/ is infinite. 

Jo 

Notice that the condition on the right-hand side is only equivalent to the non-atomicity of 
Zi - see Theorem 27.4 in [IT]. With a view towards the methods developed later in the 
present paper, let us give a short proof of the reverse inclusion in (11.11) . the direct one being 
straightforward. If u is infinite and Ti,T2 denote the two first jumping times of Z into 
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[—ri,riY, then P[T^ > 1] ^ as r] 0. Besides on {T2 < 1} one can write 
1 .1 

Jo 

and the right-hand side is a.c. on {Tg < 1} for every 77 > 0, since conditionnally on JF'' which 
is the (T-field generated by all the information given by Z on [0, 1] except T^, the variable 
{T2 — T^) has uniform hence absolutely continuous law on [0,T^], the variable AZrpv is 
non-zero a.s. and jF^-measurable, and the remaining terms are jF'^-measurable. 

This conditioning method together with, roughly speaking, a derivation procedure along 
certain jumping times, had been systematically developed in the monograph [4J where various 
absolute continuity results and smoother properties were established for several functionals of 
Wiener and other Levy processes, such as Lp-norms, proper integrals along a given function, 
one-sided and two-sided suprema. In P, 17!, it was also applied to a class of real stochastic 
equations with non-linear drift driven by Z. 

In this paper, we will deal with Non-Gaussian multidimensional Ornstein-Uhlenbeck pro- 
cesses, which are solutions to the S.D.E. 

(1.2) dXt = AXtdt + BdZt 

where A is a real n x n matrix, B a real n x d matrix and Z a rf— dimensional Levy process 
without Gaussian part. Adaptating without difficulty the discussion made in [11] pp. 104- 
105 to the case where A is not necessarily a multiple of the identity matrix, we see that the 
solution to fll.2p is given in terms of some Levy integral: 

(1.3) Xt = e'^x + [ e^'-'^^BdZs, t > 0, 

where x G is the initial condition. Apart from a natural extension of the Langevin 
equation, these processes have their roots in distribution theory because in the ergodic case 
their limit laws are known [13J to be operator-self decomposable - see [12] for further results. 
Nowadays, Non-Gaussian OU processes are also quite popular in modelling [2]. 

To state our result, we need some notation. We will assume that the reader is familiar 
with basic properties of Levy processes and jumping measures which can be found at the 
beginning of the two monographs [3],|TT]. Setting {i?^, t > 0} for the M*^- valued Levy process 
{BZt, t > 0} and z/^ for its Levy measure, we introduce the vector spaces 

Bt = Vect[ABs, s <t] and At =<A,Bt>, t > 0, 

where here and throughout, for every vector subspace S CW^ with basis {ei, . . . , Cp}, we use 
the notation 

< A,£ > = Vect[A'-^ej, i = l...n,j = l...p\. 

Notice that actually At = Vect[A^^^ABs, i = 1 . . . q, s < t], where q stands for the degree of 
the minimal polynomial of A. Setting B = ImB G M", the condition < A, B > = M" or, in 
an equivalent matrix formulation. 



(1.4) 



Rank [B,AB,...,A''-^B] = n, 
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is well-known as a controllability condition on the deterministic linear system 

x[ = Axf + But 

where {mj, t > 0} is some input function - see e.g. Chapter 1 in and the references 
therein. When (11.41) holds, we will say that {A, B) is controllable. 

Let K denote the cyclic index of A, which is the maximal dimension of its proper subspaces. 
The condition k = 1 means that A is a cyclic matrix i.e. there exists a generating vector 
6 G M" such that (6, Ah, . . . , A^~^h) forms a basis of M". This is also equivalent to g, that 
is the minimal polynomial of A is in fact its characteristic polynomial. When k > 1, it is 
the number of the invariant factors of A, viz. the unique number of subspaces Ai C M" 
such that ^1 © ... © Ak. = K", with each Ai stable by A and each A/ Ai a cyclic operator 
whose minimal polynomial ai divides ttj-i, ai being the minimal polynomial of A itself - see 
Chapter 4 in [6] or Section 0.10 in [16] for more precisions. 

Let m = Dimi3 = Rank 5. When (11.41) holds, a result of M. Heymann - see Theorem 1.2. 
in [T6] - entails that necessarily m > n. More precisely, there exist k linearly independent 
vectors bi, . . . ,1^ E B such that 

(1.5) + ... + B'' = 

with the notation B^ = Vect Abi, . . . , A'^~^bi] for i = 1 . . . k. Actually, Heymann's result 
was originally more precisely stated, connecting the subspaces B^ to the above k, invariant 
factors of A. Nevertheless we shall not need this in the sequel. Assuming (11.41) . a sequence 
(6i, . . . , br) & B of linearly independent vectors such that B^ + . . . + B^ = with the 
above notations will be called a generating sequence of with respect to {A, B), or simply 
a generating sequence when no confusion is possible. Notice that r < m. Besides, the very 
definition of k entails that necessarily r > k as well - see the proof of Theorem 1.2 in [LGj 
for details. We now come to the central definition of this work: 

Definition. With the above notations, the Levy measure is said to exhaust M" with 
respect io A if < A,B> = M" and there exists r G [n, m] and a subspace Tir d B of 
dimension r such that < A^Tij. > = and v^{l-Lr H '^^) = +oo for every hyperplane 
7i C Tir. 

This definition is related to the conditions given in [IB] for the absolute continuity of 
multivariate infinitely divisible distributions, but it is less stringent since no arcwise absolute 
continuity is required, and since may be carried by any subspace with dimension r G 
[K,m]. Here, however, the important fact is that this subspace must be chosen with respect 
to A. Introducing finally the stopping time 

T = mf{t > 0, A = K"}, 

our result reads as follows: 

Theorem. If A is non-singular, then one has 

Xi is a.c. <^=^ r = a.s. <^=^ exhausts w.r.t. A. 

Notice that the above equivalences are time- independent, so that when Xi is a.c. then Xt 
is a.c. as well for every t > 0. In other words absolute continuity is not a temporal property 
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for Non-Gaussian OU processes, unlike Levy processes without Gaussian part. Besides, when 
Xi is not a.c. then the first equivalence entails that Xi is valued with positive probability in 
some fixed affine hyperplane of M", or equivalently that a certain one-dimensional projection 
of Xi must have an atom. This, again, contrasts with Levy processes since Zi may be 
non-atomic and not absolutely continuous - see Theorem 27.19 in [TT] . 

We stress that the variable Xi itself is infinitely divisible, i.e. it is the distribution at time 
1 of some Levy process {Yt, t > 0} valued in M". Indeed, a straightfoward extension of 
Lemma 17.1 in [11] shows that Xi is ID without Gaussian part and Levy measure 



Hence, our result yields an optimal criterion of absolute continuity for a certain subclass of 
multivariate Non-Gaussian ID distributions, which we may call the OU class. To this end, 
one can check that if satisfies any condition given in [IB], then i/^ exhausts M" w.r.t. A, 
but that the converse inclusion is not necessarily true since u'^ may not be arcwise absolutely 
continous when exhausts M" w.r.t. A. 

As we mentioned before, the variables Xt converge in law when t ^ oo to some operator 
self-decomposable or OL distribution in M" under an ergodicity assumption, that is when 
the eigenvalues of A have all negative real parts and is log-integrable at infinity [T3]. If 
in addition exhausts M" w.r.t. A, then it is easy to see that the limit distribution is 
also genuinely n-dimensional. Let us notice that the absolute continuity of non-degenerated 
OL distributions had been established in |17j. This is probably related to our result, even 
though absolute continuity and non absolute continuity properties are barely stable under 
weak convergence. To make a true connection, one would need a stronger type of convergence 
such as convergence in total variation [1], but no such result seems available in the literature. 

It follows from our definition that when exhausts M", then {A, B) is necessarily con- 
trollable. On the other hand, when (A, B) is not controllable then r = +oo a.s. so that 
from our result Xi is not a.c. In a recent paper [10] which was the starting point of this 
work, it was proved that Xi is absolutely continuous as soon as {A, B) is controllable and 
the jumping measure z/ of Z is absolutely continuous in an open neighbourhood of 0. These 
conditions entail of course that exhausts M", but they are highly non-equivalent: for ex- 
ample when A is cyclic and non-singular, then our result entails that Z might be genuinely 
one-dimensional with an infinite, and possibly discrete, jumping measure carried by some 
line in B~^{h) where 6 is a generating vector of A, nevertheless Xi will be a.c. in M". Our 
method is also very different from [10], which hinges upon a certain derivation procedure, 
made possible by the absolute continuity condition on z/, along the jumping sizes of Z. Here, 
as shortly suggested above, we will differentiate along suitably chosen jumping times, the 
price to pay being the non-singularity assumption on A. Our time-derivation procedure is 
close to the one developed in [7], whose Theorem 1.1. actually entails that Xi is a.c. when 
A is non-singular and u^lT-C^) = oo for every hyperplane Ti. C W^. But again this latter 
assumption is more stringent than our exhaustion property, as well as, to the best of our 
knowledge, all conditions given in the literature on Malliavin's calculus for jump processes - 
see the references given in [9l [71 [10] . 
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In the third section of this paper we will briefly describe what happens when the driving 
process Z has some Gaussian component. By independence and by the linearity of the 
equation (11.21) . we get an analogous result which only requires a small modification of the 
proof in the Non-Gaussian case. Then we will discuss a few examples in order to provide 
more geometrical insight on the exhaustion property. In particular we will give a complete 
description in the case n = 2. Finally, we will mention some counterexamples and open 
questions. Before this we now proceed to the 



2. Proof of the theorem 
2.1. Proof that is a.c. ^ r = a.s. Setting 

So = n^* 

which is a deterministic subspace of B by the 0-1 law, it follows from the definition of Bt 
and that of a jumping measure that necessarily z/(i3g) < +oo. In particular, P[T > 1] > 
where T = inf{t > 0, Ai?^ G Sq}. Endow B with a canonical Euclidean structure and set 
Bq for the orthogonal supplementary of Bq in i3, which may be reduced to {0} if i3o = B. 
Decomposing the Levy process {-Bj, t > 0} along the orthogonal sum B = Bq® B^ : 

Bt = + B^, t > 0, 

notice that the Levy process {B^, t > 0} is either the zero process or a compound Poisson 
process with drift coefficient, say, b^. Hence, on {T > 1}, we deduce from (11.31) that 

Xi = e\ + /'e(i-^)^di?° - /'e(i-^)Vds 
Jo Jo 

where we use the notation 6-*- = if B^ = 0. Last, writing for every t G M 
(2.1) e*^ = j2Mt)A^-\ 

k=l 

where the ipr^s are certain real functions whose exact expression is given e.g. in [6] Chapter 
5, we see that 

[ e(^-^)^ dfi° e<A,Bo> a.s. 
Jo 

Putting everything together entails that if Xi is a.c. then necessarily < A, Bq > = M". But 
from the definitions of Bq and r, this yields r = a.s. 

□ 



Remark 1. The above proof shows also that if P[r > 0] > 0, then Xi is valued in some fixed 
affine hyperplane with probability P[T > 1]. Notice that if Bq = B, then this probability is 
1 and {A, B) is not controllable, which entails that actually r = oo a.s. 
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2.2. Proof that r = a.s. =^ exhausts M" w.r.t. A. We may assume < A,E > = 
since otherwise r = oo a.s. Suppose now that does not exhaust R". By the above 
assumption there exists a hyperplane Hi C B such that u^{l-i\) < oo. If < A, 7-^i > 7^ M", 
then 

r > inf{t > 0, /\Bt G T^J} > a.s. 

If < A, 7-^1 > = M", then there exists a hyperplane I-L2 C TYi such that u^iTii fl Tig) < 00 
because does not exhaust M", and in particular 1/^(7-^2) < 00 since we already have 
ly^iWi) < 00. Similarly, when < A,7^2 > 7^ K", then 

r > inf{t > 0, /\Bt G 7^^} > a.s. 

When < y4, 1-L2 > = M", one can then repeat a finite number of times the same discussion as 
above: alltogether this entails that r > a.s. if does not exhaust M", which completes 
the proof by contraposition. 

□ 

2.3. Proof that exhausts M" w.r.t. A ^ Xi is a.c. This is the difficult inclusion 
and we will first establish three lemmas. The first one is an easy application of the implicit 
function theorem. The second one is an a.s. independence result on a certain class of linear 
systems, for which we could not find any reference in the literature on control theory. The 
third one allows to choose suitably the jumping times of B which will be later targeted 
into Xi via a certain a.s. submersion. Throughout, A will stand for the Lebesgue measure 
independently of the underlying Euclidean space. 

Lemma 2. Let X he an absolutely continuous random variable in W with p > n and 
(f : MP ^ R" be a function such that its Jacobian matrix dip verifies Rank d(f{X) = n 
a.s. Then <f{X) is absolutely continuous in M". 

Proof: Choose any Xj^j- G fl A/"^, where Af = {x & W / Rank d(f{x) < n}. Setting 

X = Xl^xeAf^} + a;_vl{XGAr}, 

we see by assumption that X = X a.s. so that ip{X) = (p{X) a.s. as well. Hence, it suffices 
to show that ip{X) is absolutely continuous in R". In the case p = n this is a well-known 
fact for which we found a proof in [H], Lemma IV.3.L For the sake of completeness, we will 
give an argument in the general case p > n. 

Fix an underlying probability space P). By approximation, it is enough to show 

that for every relatively compact set F C $ = {ip{X{u!)), u & Q} such that A(F) = 0, 

P[V9(X) G F] = 0. 

For every y E T, fix x E (p~^{y) C A. Since R.a.Bkdip{x) = n, by the implicit function 
theorem there exist Vx and Wy open neighbourhoods of x and y respectively. Op and On open 
neighbourhoods of respectively in MP and R" endowed with a canonical basis, ipx '■ ^ Op 
and ipy : Wy On diffeomorphisms such that 

ipy o ip o : Op ^ On 
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is the canonical projection from Op to On- Taking a finite covering {W^^, . . . , Wj^^,} of V 
yields 

k 



i=l 
k 



i=l 



with the notation Fj = ipyX^yi l~l T), which has Lebesgue measure zero in since ipy. is 
a diffeomorphism. Setting Ai = {'ipy. o (p{X) G Fj} for every i G {!,..., k}, the random 
variables Yi = Ij^^-ip^-lX) are absolutely continuous in since tp^i are diffeomorphisms. 
Besides, since the projections ipy. oipoil)^^ have full rank, from Lemma 2.3 in the random 
variables ipy^ o {p o ipxiiyi) are absolutely continuous in . Hence, for every i G {1, . . . , k}, 

F[^y^ o G F,] = p [ijy^ o ^ o ij.xyi) e r,] = 0, 

which entails that P[(y9(X) G F] = and completes the proof. 

□ 

Lemma 3. Assuming ( [i.^p , let [bi, . . . ,br) be a generating sequence with respect to {A,B) 
for some r G [K,m]. Then the set 



|(t},...,t;,...,t;,...,gGM^^VRank 

has zero Lebesgue measure. 



e*i^6i, . . . , e*9^6i, . . . , e^^^\, . . . , e*^^6. 



< n 



Proof: We first consider the case k = r = 1. Fix 6 G a generating vector such that 
= Vect [b, Ab,..., A'^-^b] . The function 

^ Det [e''\...,e'-H] 

is analytic in M" and it is not identically zero. Actually, if it were, then the analytic function 
p : t ^ e^^b would be valued in some fixed hyperplane of M", as well as all its successive 
derivatives, which is impossible since 

Vect [p(0),p'(0),...,p("-^)(0)] = Vect [b, Ab, A'^-'b] = M". 

Hence, since the zero set of a real analytic function over M" either is itself or has zero 
Lebesgue measure - see [S] p. 240 or Lemma 2 in [17], we obtain that 



Rank [e*i^6,...,e*"^6] 



n 



almost everywhere in M", which completes the proof when k, = r = 1. 

We now proceed to the remaining cases. Recalling (12.11) . we first claim that the q x q 
matrix 
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has rank q a.e. in R'^. Indeed, by the definition of q there exists 6 G such that Rank 
[b, Ab, . . . , A'^~^b] = q. Setting Ab = Vect[6, Ab, . . . , A'^^^b] and viewing A as a cychc endo- 
morphism of the g— dimensional vector space At, we see from the case k = r = 1 that Rank 
[e^'^^b, . . . ,e'^''^b] = q a.e. in W. However, it follows from (12.11) that 



e 



*^^6,...,e*'^6] = [b,...,A'^-'b]x^iigit,,...,tg 



SO that ^q{ti, . . . ,tq) must have rank q a.e. in W as well. Let now . . . , br) be a generating 
sequence with respect to {A, B). Setting Bi = Vect[6i, Abi, . . . , A'^~^bi], we have Dim Bi < q 
for every i = 1 . . .r. Similarly as above. 



J^Al taA 



-'11 ■ ■ ■ 1 



-X] = [6„...,A^-i6,]xvl/^(ti,...,g 



and since \E'q is a.e. invertible, this entails that Rank [e'^^^foj, . . . ,e*'''^6j] = Dim Bi a.e. in 
W. Notice that by fl2.ip . one has e*^6j G Bi for every t G M, so that [e^^^foj, . . . , e^'^foj] forms 
actually a basis of Bi a.e. in M'', for every i = 1 . . .r. Besides, it follows from the definition 
of the Lebesgue measure that if ^i, . . . ,Ar are negligible sets in M'^, then [Al x . . . x Af)'^ 
is negligible in Mfl^"^ . Putting everything together entails that a.e. in M'^^'', 



Rank 



3*i%,...,e*J%,...,e*^X...,e*SX 



Dim(-Bi + . . . + S^) = n, 



where the last equality comes from the definition of the generating sequence . . . , fe,.)- The 
proof is complete. 

□ 

Remarks 4. (a) By the same argument, one can prove that 

A { (ti, . . . , t„) G / Rank [e^^^E, . . . , e*"^5] < = 

as soon as fll.4p holds. An interesting point is that when A has real spectrum, then fll.4l) 
actually entails that 

(2.2) Rank [e*i^5, . . . , e*"^5] = n 

for every ti < ... < t„. The latter is false nevertheless when A has non-real eigenvalues. 
The proofs of the above two facts involve rather technical considerations on generalized 
Vandermonde matrices, which we shall not discuss here. 

(b) I could not find in the literature any material on the following question, whose positive 
answer would quickly entail (12.21) . Assuming that A has real spectrum and that (II. 4p holds, 
let u^, . . . , be n piecewise constant control functions which are linearly independent over 
[0, 1] and consider the linear systems 

= X + j {Axl + Bui) ds, t > 0, i = 1 . . . n. 

Does {x\, . . . , Xi) form then necessarily a basis of M*^? 

A family (Ci, . . . , Cr) of disjoint pointed cones with common vertex at zero such that every 
(ci, . . . ,Cn) G (Ci, . . . , Cr) is a generating sequence with respect to {A, B) will be called a 
generating garland with respect to {A, B), or simply a generating garland when there is no 
ambiguity. When (6i, . . . ,br) is a generating sequence, notice that (Ci, . . . , Cr) defined by 



ABSOLUTE CONTINUITY OF ORNSTEIN-UHLENBECK PROCESS 



9 



Ci = {fibi, /i > 0} for i = 1 . . . r, is a generating garland. Our last lemma makes a connection 
between this notion and the exhaustion property: 

Lemma 5. If i'^ exhausts M" w.r.t. A, then for every M > there exists r G [k, m] and a 
generating garland (Ci, . . . , Cr) such that v^{Ci) > M for every i = 1 . . .r. 

Proof Fix M > 0. If exhausts R" w.r.t. A, then we know that < A,B >= M". 
Suppose first that viT-f^) = oo for every hyperplane TC C B and let us show that there 
exists a generating garland (Ci, . . . , Cm) such that v^{Ci) > M for every i = 1 . . .m, which 
is intuitively obvious. 

If denotes the unit Euclidean sphere of B, consider a family {II^, 5 > 0} of finite 

measurable partitions of S"^'^ such that Diam(n5) < 6 and Us' is a subpartition of II^ for 
every 5' < 5. Let denote the disjoint finite family of pointed cones with vertex at zero 
and apex in 11^ such that v^iC) > M for every C G C^^. If for every 5 > no generating 
garland of size m is contained in C£f , then for every 5 > there exists at least one hyperplane 
Tis intersecting every C G Cf^, and the assumption Diam(n5) < 5 readily entails that all 
these 7^5 's converge - in the sense that their normal unit vectors converge in the metric 
space S'^~^ - to some fixed hyperplane TIq. Last, it is a bit tedious but not difficult to see 
that by construction and by the finiteness of II^, one must have v{T^q) < oo, which yields a 
contradiction. Hence, there exists 6 > such that C|,^ contains a generating garland of size 
m, and we are done. 

If viT-C^) < oo for some hyperplane Ti. G B, then by definition of the exhaustion property 
there must exist an hyperplane Hm-i such that < A,T-Cm-i > = R"- If ^{'Hm-i H H'^) = oo 
for every subspace Ti. C Tim-i-, then reasoning exactly as above we can show that there 
exists a generating garland (Ci, . . . , Cm-i) such that v^iCi) > M for every z = 1 . . . m — L If 
^{'Hm-i nTi'^) < oo for some subspace H C Hm-i, then we can repeat the same procedure as 
above. But again by the definition of the exhaustion property, the latter procedure cannot 
be repeated more than (m — n) times: all in all, this shows that there exists r G [k, m] and 
a generating garland (Ci, . . . , C,.) such that ^^{Ci) > M for every i = 1 . . . r. 

□ 

Remark 6. When exhausts M", there may exist no generating garland (Ci, . . . , C^) such 
that v^{Ci) = oo for every i = 1 . . .r. Think of the situation where k = m = 2 and is 
infinite with support in the arc y = x'^, B being endowed with an orthonomal frame Oxy. 

End of the proof. From f 1 1.3 1) and after time-reversal, it is enough to prove that 

Y = [ e'^BdZ, = [ e^^dE, 
Jo Jo 

is absolutely continuous. For this purpose, we will use the same method as depicted in the 
introduction, in a somewhat more elaborated manner. If F C M" is such that A(F) = 0, we 
need to show that for every e > 

P[FGr] < e. 

Fix e > and let M > be such that 



¥[T%, > 1] < s/m 
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where T^li is the sum of (g + 1) independent exponential variables with parameter M. By 
LemmaO there exists r G [n, m] and a generating garland (Ci, . . . , C,.) such that v^{Ci) > 2M 
for every i = 1 . . .r. Besides, if B^j stands for the Euclidean ball of B centered at with radius 
rj and if = Ci (1 B"^, then we can actually choose rj > such that v^iC^) > M for every 
i = 1 . . .r. Let {T^^p, p > 1} be the ordered sequence of jumping times of the Levy process 
{-Bf , t > 0} into and set = sup{T^^p, i = 1 . . . r} for every p > I. We have 

r 
i=l 

and it is hence sufficient to prove that 

(2.3) p [t;+i < 1, F G r] =0 

for every r] > 0. Let JF'? be the a— algebra generated by {T^'p, p > q + I, i = l...r}, 
{ABrpv , p > 1, i = 1 . . .r} and the Levy process defined by 

Bl = Bt-^ ABtv , t>0. 
On {Tqj^i < 1}, one can write 

Y = Y^ + ^^e^-^Ai?^.^. 
i=i j=i 

with a JF^— measurable random variable. Since Tj'^^ is JF*?— measurable as well, we have 

p [r;Vi < 1, r G r] = p [r;Vi < i, p [r g r | j^*?]] 

= P[T;'+i < 1, P[f'' G F'' I r'']] 
where = F — y is a JF''— measurable set such that A(F'') = 0, and 

= ^^e^-^AE^.^.. 

i=i j=i 

The key-point is that by standard properties of jumping measures, since the C^^'s are disjoint, 
conditionally on JF*? the law of the (r x g)— uple (T^^, . . . , T^^, . . . , Tjjj^, . . . , T^^) is that of 
the tensor product of r independent g-th order statistics respectively on [0,T^^g_^_l], viz. the 
tensor product of r independent uniform laws on the respective sets 

{o<n<...<t;<7::,+i}. 

In particular, the law of (T^'^i, . . . , T^^^, ■ ■ ■ , T^i, ■ ■ ■ , T^g) is absolutely continuous in and 
by Lemma [2l (12. 3p will hold as soon as the Jacobian matrix of the application 



1=1 j=i 
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from M'^^'' to M" has rank n a.e. conditionnally on JF*? (recall that the sequence {^Brpn , p > 

1, i = 1 . . . r} is jF^-measurable and independent of {^i^p, P < 9, i = 1 ■ ■ ■ r}). This Jacobian 
matrix is equal to 

(2.4) A X pi^AS^.^, . . . ,e*J^A5r.^, . . . , e*^^A5^.^ , . . . ,e*S^Ai?r.^ 

and, since A is invertible and AZS,^ G for every i = 1 . . . r and j = 1 . . . g, conditionnally 
on J-"^ it has a.s. the same rank as 

e*i^6i, . . . , e*9^6i, . . . , e^i^fe^, • • • , e^'^&r 

where (6i, . . . , 6^) is some generating sequence of M" with respect to (A, 5). Now by Lemma 
131 the latter has full rank a.e. and the proof is finished. 

□ 

Remark 7. The invertibility assumption on A is only useful for getting the full rank a.s. 
of the Jacobian matrix given in (12. 4p . Nevertheless this is a crucial assumption, and in the 
next section we will give a counterexample when A is singular. 

3. Final remarks 

3.1. The case with a Brownian component. If the driving Levy process Z has a non- 
trivial Gaussian component, then by the linearity of (11. 2p this amounts to consider the 
problem of absolute continuity for 

Xt = e'^x + [ e^'-'^^dWs + [ e(*-^)^d5„ t>0, 
Jo Jo 

where W is some i3-valued Brownian motion independent of the Non-Gaussian Levy process 
B. Set H = < A,lmW > and, given any Euclidean structure on M", denote by H-^ the 
orthogonal complement of if in M". The i7— valued random variable 

"1 







is Gaussian and by a classical result in control theory - see e.g. Theorem 1.1 in (TB] - it 
is non-degenerated, hence absolutely continuous in H. Since W and B are independent. 
Lemma 3 in [17] and Lemma 2.3 in [10] yield 

Xiisa.c. ^ -^^^ (/ is a.e. 

where Hh^ stands for the orthogonal projection operator onto H-^. A straightforward mod- 
ification of our proof entails then 

Xiisa.c. <^==^ Tjj = a.s. <^=^ z/'^ exhausts w.r.t. {A,H). 

where with the notations of the introduction we set t^j = infjt > 0, At + H = M."'}, n^j for 
the minimal number of linearly independent vectors bi, . . . ,bp & B such that 

< A,bi> +...+ < A,bp> + = M", 
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and say that exhausts with respect to (A, H) ii < A,B > + if = and there exists 
r G [n^^m] and a subspace Tir G B of dimension r such that < A^Ti^ > + H = W"' and 
u^iHr n T-C^) = +00 for every hyperplane H C Hr- 

3.2. Some explicit descriptions of the exhaustion property. From now on we will 

assume that Z has no Gaussian part, that is B has no Gaussian part either. Let us first 
consider the case n = 1, i.e. X is solution to 

(3.1) dXt = aXtdt + dBt 

where a G M and B is one-dimensional. The exhaustion property just means that is 
infinite and our result reads 

(3.2) Xi is a.c. <^=^ is infinite 

as soon as a 7^ 0. Notice that this is actually an immediate consequence of Theorem A in 
[0] - see also Theorem 1.1. in [7j. Let us also give a short proof of the non-trivial reverse 
inclusion in (13.21) . similar to that given in the introduction: the solution to fl3.ll) is given by 

Xt = e*"x + [ e'^^*"") dBs = e*"x + Bt + a [ e<^-''^B, ds, t > 0, 
Jo Jo 

where x is the initial condition and where in the second equality we made an integration by 
parts, assuming _Bo = without lost of generality. Hence, leaving the details to the reader, 
one may follow roughly the same method as for the integral of B, noticing with the same no- 
tations that on {T2 < 1} the value of Bi does not depend on T^, hence Bi is jF^'-measurable 
as well. 



Let us now discuss the case tt, = 2. To simplify the notations, we will denote by X the set 
or family of sets where is infinite if and only if the exhaustion property holds. We will 
also suppose implicitly that A is non-singular. Up to some equivalent transformations on A 
which are not relevant to the absolute continuity problem, there are four situations: 

(a) A has no real eigenvalue, in other words A is a multiple of 

/ cos 6' sin^^ \ 
\^ — sin 9 cos 9 J 

for some 9 g]0, 7i[. Then k = 1 i.e. A is cyclic, and it is easy to see that every non-zero vector 
of is generating. Hence we simply have X = viz. as in the real case, Xi is a.c. if and 
only if z/^ is infinite. 

(b) A is a multiple of the identity matrix. Then k, = 2 and < A,b > = Vect{6} for every 
b G M^. This means that X = {(Vect{6})'^, b G M^} viz. the infinite part of must not be 
carried by any line in M^. 

(c) A is a Jordan cell matrix, i.e. A is of the type 

(a l\ 
\0 a J 
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with q; 7^ 0. Then k — 1 and every non-zero vector of M? is generating except the multiples 
of (1, 0) : we have I = (Vect{(l, 0})^ 

(d) A = Diag(a,/3) with a ^ f3 and a,f3 non zero. Then k, = 1 and every non-zero vector 
of is generating except those in Vcct {(1,0)} U Vcct{(0, 1)}. On the other hand, (Vect 
{(1, 0)} - {0}, Vect {(0, 1)} - {0}) is a generating garland: we have J = {(Vect{(l, 0)})^ n 
(Vect{(0, 1)})'^} U {Vect{(l, 0)} and Vect{(0, 1)}}. 

When n > 2, it becomes quite lengthy to depict the exhaustion property. Let us give 
however four typical examples when n = 3, keeping for X the same meaning as above and 
using the notations Hx = {x = 0}, Tiy = {y = 0} and Ti^ = = 0} where Oxyz is a given 
orthogonal frame of M?. 

(f ) A is a Jordan cell matrix, i.e. A is of the type 



with a 7^ 0. Then k — 1 and every non-zero vector of is generating except those in H.^ : 
we have X = 7i^. 

(g) ^4 is a block matrix of the following type 



with a ^ (3 and a,/? non zero. Then k — 1 and every non-zero vector of is generating 
except those in Hx^Hz '■ we have X — 'W^n'W^. 
(h) A is a block matrix of the following type 



with a 7^ 0. Then k = 2 and every generating sequence must not be valued in any hyperplane 
Tiu = W"*- with u a unit vector of Oxz : we have X = {ht^, u G Oxz}. 

(i) A = Diag(Q;, /3, 7) with distinct non zero a, {3 and 7. Then n = 1 and every vector in 
'W^m-LyH HI is generating. But as in dimension 2, one can also build generating garlands 
with one component in Tix^'Hy or Tiz- The infinity set X is then the union of the following 



eight sets or families of sets: HlnHinHl, {HxnninHl andUynHl}, {T^j^nT^^nH^ andT^^n 

ni},{nynn'',nn'x andTi^ n7f^},{H, n^;; n^^ andH^ n7i^},{H, n^^ n?i^ andT^^ n 
W;}, {Hx nninm andH^ n ni], {Hx nHyDH'; andUy n n andH^ n 7^^: n m}. 



3.3. Some open questions. As we mentioned before, our theorem no more holds when A 
is singular, as shows the following counterexample with n = 2,d — m — 1, 







and 
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From the control theory viewpoint, this example yields the so-called rocket car equations, 
which serve as toy-models [S] for studying the Pontryagin maximum principle. From the 
stochastic viewpoint, assuming x = (0,0) one gets the so-called Kolmogorov process X = 
{X\X^), with 



Xl = Zt and X^ = Z^ds, t> 0, 

Jo 

and Z is a one- dimensional Levy process. Notice that {A, B) is controllable, so that 
exhausts w.r.t. A if and only if is infinite. But then X\ = Z\ may be singular - see 
again Theorem 27.19 in |11J - so that by Lemma 2.3 in [TD], X^ is not absolutely continuous 
either. 

When A is singular and m, one may wonder if the following holds true for every t > : 

Bt is a.c. in R" =^ Xt is a.c. in M". 

From our theorem, this property is trivial when n = 1 and we also refer to Theorem B 
in ^ for a non-linear extension. When n > 2 some geometrical difficulties arise however, 
which will be the matter of further research. When m < n the problem seems much more 
complicated without further conditions on the jumping measure. In particular I do not have 
the answer to the following basic question, which would solve the problem at least for the 
Kolmogorov process: 

For a real Levy process {Zj, t > 0}, if Zi is a.c. in R, is ^Zi, J Zfdt^ a.c. in R^? 

To conclude this paper, let us go back to the case n = 1 and consider the following class 
of infinitely distributions 



O = e'^ dB^ , B real Levy process j- 



for which we proposed the name OU class. If fi^ G O corresponds to some Non-Gaussian 
Levy process B, after time-reversal our previous discussion entails 

(3.3) fx^ is a.c. <^=^ is infinite. 

Besides, with a little refiexion, one can show that (13.31) also holds when replacing in O the 
kernel e* by any function f{s) whose derivative does not vanish in ]0, 1[. In particular, 
the equivalence (13. 3p will hold for the well-known U-class where f{s) = s, and B-class where 
/(s) = logs. Of course, for these two special classes one could get (13. 3p directly in considering 
the jumping measure of /i^ which happens to be absolutely continuous - see [I] for details - 
and applying Tucker's result - see Theorem 27.7 in [TT]. It would be interesting to investigate 
which exact class of integration kernels entails (13.31) for Levy integrals, and also what occurs 
in the multivariate case. 
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